Accepted by The Astronomical Journal 

Preprint typeset using IATgX style emulatcapj v. 04/03/99 



O 
O 

o 

(N 

c 

a 

(N 



> 

in . 

o : 
o . 
o ■ 

On; 

6 



13 



A SEARCH FOR SUB-MILLISECOND PULSATIONS IN UNIDENTIFIED FIRST AND NVSS 

RADIO SOURCES 

Fronefield Crawford 1 and Victoria M. Kaspi 2 

Department of Physics and Center for Space Research, Massachusetts Institute of Technology, Cambridge, 

MA 02139, USA 
AND 

Jon F. Bell 3 

The University of Manchester, Jodrell Bank, Macclesfield, Cheshire, SK11 9DL, UK 
Accepted by The Astronomical Journal 

ABSTRACT 

We have searched 92 unidentified sources from the FIRST and NVSS 1400 MHz radio survey catalogs 
for radio pulsations at 610 MHz. The selected radio sources are bright, have no identification with 
extragalactic objects, are point-like and are more than 5% linearly polarized. Our search was sensitive 
to sub-millisecond pulsations from pulsars with dispersion measures (DMs) less than ~ 500 pc cm -3 
in the absence of scattering. We have detected no pulsations from these sources and consider possible 
effects which might prevent detection. We conclude that as a population, these sources are unlikely to 
be pulsars. 

Subject headings: pulsars: general — surveys — radio continuum 



1. INTRODUCTION 

The FIRST survey (Faint Images of the Radio Sky at 
Twenty Centimeters) and NVSS survey (NRAO VLA Sky 
Survey) are recent 1400 MHz VLA radio surveys of the 
Northern sky. The FIRST survey is an ongoing survey of 
the North and South Galactic caps using the VLA in B- 
configuration with a synthesized beam size of 5.4" (Becker 
et al. 1995). In the published FIRST catalog of radio 
sources from the first two observing sessions in 1993 and 
1994 (White et al. 1997), 1550 square degrees of the North 
Galactic cap were covered spanning 7 h < a(J2000)< 18 h 
and +28° < <5(J2000) < +42°. The positions and flux 
densities of ~ 1.4 x 10 5 discrete radio sources 4 are com- 
plete down to a flux density of ~ 1 mJy. The NVSS survey 
(Condon et al. 1998) covers 5 > -40° (covering 82% of the 
celestial sphere) and catalogs more than 1.8 x 10 6 sources 
complete down to a flux density of ~ 2.5 mJy. The NVSS 
survey was conducted with the VLA in D and DnC config- 
urations with a synthesized beam size of 45". The NVSS 
survey also preserves polarization information. 

Several large-scale pulsar surveys have previously been 
conducted at high Galactic latitudes (see Camilo 1997 and 
references therein). However, the rates at which the re- 
ceived analog power was sampled and digitized in these 
surveys, typically 3-4 kHz, and the low observing radio 
frequencies (~ 400 MHz), combined with relatively large 
radio frequency channel bandwidths of between 125 kHz 
and 1 MHz, restricted their sensitivity to sub-millisecond 
pulsars to very small dispersion measures (DMs) (DM < 
10 pc cm -3 ). Large-scale surveys that maintain sensitiv- 
ity to sub-millisecond periodicities over a wide range of 
DMs are difficult: the fast sampling rate and small radio 



frequency channel bandwidth required make large total 
bandwidths and long integration times currently imprac- 
tical. However, a targeted search for sub-millisecond pul- 
sations is possible using narrow frequency channels and a 
fast sampling rate. Such a survey is of course also sensi- 
tive to long-period pulsars which may have been missed 
in previous surveys due to radio frequency interference or 
scintillation. 

Consideration of the properties of known recycled pul- 
sars and representative models of magnetic field decay and 
equations of state suggests that a significant population 
of sub-millisecond pulsars could be present in the Galaxy 
(e.g., Possenti et al. 1998). It is possible, therefore, that 
some of the sources which remain unidentified in radio sur- 
vey catalogs could be bright sub-millisecond radio pulsars 
which have previously escaped detection in high-latitude 
pulsar surveys. To date, no pulsar has been found with 
a period shorter than that of the first millisecond pulsar 
discovered, PSR B1937+21, which has a 1.56 ms period 
(Backer et al. 1982). The discovery of a sub-millisecond 
pulsar would place important constraints on the equation 
of state of neutron matter at high densities (e.g., Kulkarni 
1992). 

2. TARGET CHOICE AND OBSERVATIONS 

The FIRST and NVSS surveys contain a number of 
bright sources which are unresolved and have no identi- 
fication in other wavebands. Although over 99% of bright 
sources (S1400 > 60 mJy) found in previous large-scale sur- 
veys are believed to be active galactic nuclei (AGN) (Con- 
don et al. 1998), many sources in the FIRST and NVSS 
catalogs remain unidentified. One possibility is that they 
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are previously unrecognized radio pulsars. Since pulsars 
often have a high degree of linear polarization (Lyne & 
Manchester 1988), polarized sources are good targets for 
pulsar searches. Han & Tian (1999) have identified 97 ob- 
jects in the NVSS catalog which are coincident with known 
pulsars. Of the 89 redetected pulsars in Table 1 of their 
paper for which the degree of linear polarization could be 
determined from the NVSS observations, only 8 had an ob- 
served nominal fractional linear polarization less than 5%. 
The intrinsic degree of polarization of these pulsars may 
even be higher if bandwidth depolarization effects are sig- 
nificant, in which case an even higher fraction of the pulsar 
sample is more than 5% linearly polarized. Figure 2 of Han 
& Tian (1999) shows that only ~ 10% of identified quasars 
and ~ 10% of BL-Lac objects are more than 5% linearly 
polarized. Thus, although there is not a clear polarization 
cutoff separating the pulsar and extragalactic populations, 
a polarization threshold of 5% excludes most (~ 90%) of 
the identified non-pulsar population while retaining the 
majority (~ 90%) of the identified pulsar population. 

We have searched for radio pulsations in bright (6*1400 > 
15 mjy) point-like unidentified sources from the FIRST 
and NVSS survey catalogs which are more than 5% lin- 
early polarized at 1400 MHz. Sources were selected di- 
rectly from the catalogs if they met certain criteria. 

Unidentified FIRST sources had their flux densities 
checked against their corresponding NVSS flux densities. 
If a source were extended, the better resolution of the 
FIRST survey would be expected to yield a lower flux 
density for the source than the NVSS survey. Therefore, 
in order to eliminate extended objects, sources were only 
included if their FIRST and NVSS flux densities agreed 
to within a few percent (indicating an unresolved non- 
variable source) or if the FIRST flux density exceeded the 
NVSS flux density (indicating an unresolved scintillating 
source) . The large number of extended sources in the cat- 
alogs makes this filter necessary, though it does unfortu- 
nately eliminate scintillating sources which happen to be 
fainter in the FIRST survey. 

For unresolved NVSS sources outside of the FIRST sur- 
vey region, pointed VLA observations were undertaken in 
October 1995 in B-configuration in order to obtain the 
same angular resolution (5.4") as the FIRST survey (R. 
Becker & D. Helfand, unpublished work). Sources in these 
observations were then subjected to the selection criteria 
described above. A total of 92 objects from the catalogs 
(39 appearing in both FIRST and NVSS, and 53 appear- 
ing only in NVSS) fit our selection criteria and were suffi- 
ciently far north to be observed in our search. The posi- 
tions, NVSS total intensity peak flux densities, and frac- 
tional linear polarization from the NVSS peak flux values 
are listed for the selected sources in Tabic 0. 

Each of the 92 unidentified sources was observed at a 
center frequency of 610 MHz in two orthogonal linear po- 
larizations for 420 s with the Lovell 76-meter telescope 
at Jodrell Bank, UK. A total bandwidth of 1 MHz was 
split into 32 frequency channels with detected signals from 
each channel added in polarization pairs and recorded on 
Exabyte tape as a continuous 1-bit digitized time series 
sampled at 50 /is. 

The minimum detectable flux density of periodicities in 
a pulsar search depends upon the raw sensitivity of the 
system and a number of propagation and instrumental ef- 



fects (e.g., Dewey et al. 1985). Interstellar dispersion con- 
tributes to the broadening of the intrinsic pulse according 
to 



TDM = 




Here t^m is in milliseconds, v is the observing frequency 
in MHz, DM is the dispersion measure in pc cm~ 3 , and 
A.1/ is the channel bandwidth in MHz. For our observing 
system, tdm is 1.135 /is per pc cm~ 3 of DM. The fast sam- 
pling rate (t S amp = 50 /its) and small channel bandwidth 
(Az/ = 31.25 kHz) in our survey made it sensitive to sub- 
millisecond pulsars for a large range of DMs (DM < 500 
pc cm -3 ) in the absence of pulse scattering effects. Our 
estimated sensitivity to pulsations for a range of periods 
and DMs is shown in Figure |l|. A detailed explanation of 
the calculation used to produce Figure [I] can be found in 
Crawford (2000). 

3. DATA REDUCTION 

In each observation, the frequency channels were dedis- 
persed at 91 trial DMs which ranged from to 1400 pc 
cm' 3 and were summed. We searched this large DM range 
in the unlikely event that a source could be a previously 
missed long-period pulsar with a high DM. For DM > 500 
pc cm -3 the channel dispersion smearing is too great to 
maintain sensitivity to sub-millisecond pulsations from all 
our sources. However, the Taylor & Cordes (1993) model 
of the Galactic free electron distribution indicates that for 
all of our sources, with the exception of two that are within 
5° of the Galactic plane, the DMs are expected to be less 
than 100 pc cm' 3 regardless of distance. Each resulting 
dedispersed time series of 2 23 samples was then coherently 
Fourier transformed to produce an amplitude modulation 
spectrum corresponding to a trial DM. 

Radio frequency interference (RFI) produced many false 
peaks in certain narrow regions of the modulation spectra 
at low DMs. We therefore masked several frequency ranges 
and their harmonics in which RFI appeared regularly so 
that any true pulsar signals were not swamped by interfer- 
ence. Typically 1-2% of the modulation spectrum in each 
observation was lost in this way. 

We then looked for the strongest peaks in the spec- 
tra. First, each modulation spectrum was harmonically 
summed. In this process, integer harmonics in the modu- 
lation spectrum are summed, enhancing sensitivity to har- 
monic signals (e.g., Nice, Fruchter, &; Taylor 1995). This 
is particularly useful for long-period pulsars, which have 
a large number of unaliased harmonics. After summing 
up to 16 harmonic signals, the highest candidate peaks 
in the modulation spectrum were recorded along with the 
period, DM, and the signal-to-noise ratio (S/N). Redun- 
dant harmonic candidates were then eliminated. Unique 
candidates were recorded if they had S /N > 7 and if the 
candidate period appeared in at least 10 DM trials. The 
final candidates for each beam were inspected by dedis- 
persing the original data at DMs near the candidate DM 
and folding the data at periods near the candidate period 
in order to look for a broad-band, continuous pulsar-like 
signal. 

This technique was tested by observing several known 
bright pulsars (PSR B1937+21, PSR B0329+54, and PSR 
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J2145— 0750) throughout the survey. The results for these 
pulsars are listed in Table ||. All three pulsars were de- 
tected with S/N consistent with our survey sensitivity, 
though scintillation affects the detection strengths. 

4. DISCUSSION 

We did not detect any significant pulsations from the 
target sources. Here we consider possible effects which 
could prevent detection if they were pulsars. 

The selected sources were bright, with the weakest 
source having a 1400 MHz flux density of 15 mJy. As- 
suming a typical pulsar spectral index of a = 1.6 (Lorimer 
et al. 1995), where a is defined according to S ~ v~ a , this 
source would have flux density 58 mJy at 600 MHz (the 
horizontal dashed line in Figure [l]). For expected DMs, 
this is about seven (four) times greater than our sensitiv- 
ity limit for periods greater than (about equal to) 1 ms, as 
indicated in Figure [|. All of our sources, therefore, were 
bright enough (in the absence of scintillation) to be easily 
detectable with our observing system if they were pulsars. 

Dispersion smearing is not a factor preventing detec- 
tion, since all but two of these sources have high Galactic 
latitudes (\b\ > 5°) and should have DM < 100 pc cm -3 
regardless of distance. This is well within our sensitivity 
limits to sub- millisecond pulsations (see Figure [l]). Inter- 
stellar scattering, which can be estimated from the Taylor 
& Cordes (1993) model of the Galactic electron distribu- 
tion, is expected to be negligible at 610 MHz. Two of the 
sources (0458+4953 and 0607+2915) are within 5° of the 
Galactic plane; for an assumed DM of 100 pc cm -3 , their 
predicted pulse scatter-broadening times are ~ 140 [is at 
610 MHz. This is small enough to maintain sensitivity to 
sub-millisecond pulsations, and is of the same order as the 
dispersion smearing (113 [is) at this DM. 

An extremely wide beam would prevent modulation of 
the pulsed signal and could render a pulsar undetectable. 
However, this would likely occur only in an aligned rota- 
tor geometry in which both the spin and magnetic axes 
are pointing toward us. This is unlikely for our targets, 
since, if they were pulsars, the position angle of linear 
polarization is expected to follow the projected direction 
of the magnetic axis as the star rotates (Lyne & Manch- 
ester 1988). This geometry would significantly reduce the 
measured degree of linear polarization as the pulsar ro- 
tates, inconsistent with our choice of significantly polar- 
ized sources. 

Scintillation is the modulation of a radio signal pass- 
ing through a medium of variable index of refraction, such 
as an inhomogeneous interstellar plasma. In diffraction 
scintillation, the wave scattering causes interference which 
can enhance or suppress the amplitude of the radio signal 
on the time scale of minutes (e.g., Manchester & Taylor 
1977). These intensity fluctuations vary as a function of 
radio frequency at any given time and have a characteristic 
bandwidth A^, where 

Azz-llz, 22 / 5 ^ 11 ^. (2) 

Here Aj/ is in MHz, v is the observing frequency in GHz, 
and d is the distance to the source in kpc (Cordes, Weis- 
berg, & Boriakoff 1985). For pulsars with d < 1 kpc, this 
characteristic bandwidth exceeds our 610 MHz observing 
bandwidth of 1 MHz. Indeed, of 28 nearby pulsars ob- 



served at 660 MHz by Johnston, Nicastro, & Koribalski 
(1998), 13 had scintillation bandwidths greater than our 
bandwidth of 1 MHz and had a characteristic fluctuation 
time-scale greater than our integration time of 420 s. More 
than half of these 13 pulsars had distances less than 1 kpc, 
and only one had d > 2 kpc. If the sources we surveyed 
were placed at a distance of 2 kpc and a spectral power 
law index of a = 1.6 is assumed, their 400 MHz luminosi- 
ties would all be at the upper end of the observed pulsar 
luminosity distribution (£400 > 450 mJy kpc 2 ). This sug- 
gests that if these sources were pulsars, they are likely to 
be closer than 2 kpc and therefore in the distance range 
where the scintillation bandwidth exceeds our observing 
bandwidth. In this case, scintillation causes the probabil- 
ity distribution of the observed intensity to be an expo- 
nential function with a maximum in the distribution at 
zero intensity (McLaughlin et al. 1999). The number of 
sources we expect to sec in our sample is the sum of the 
probabilities that we will see each individual source (i.e., 
that the scintillated flux is above the minimum detectable 
flux). For the range P > 1 ms, we expect to see 88 of 
the 92 sources (5% missed). For the range P ~ 1 ms, we 
expect to see 85 of the 92 sources (8% missed). Thus, only 
a few of our sources are likely to have been missed due 
to scintillation. Therefore scintillation cannot account for 
the non-detection of the bulk of the sources in the survey. 

A pulsar in a binary orbit experiences an acceleration 
which changes the observed modulation frequency during 
the course of the observation. Sensitivity to pulsations is 
degraded if the frequency drift exceeds a single Fourier bin 
A/ = 1/Ti n t, where T; nt is the integration time. Assuming 
that the acceleration of the pulsar is constant during the 
observation, a critical acceleration can be defined, above 
which the change in frequency is greater than A/ and the 
sensitivity is reduced: 

Ocrit = 77p2 ■ (^) 
J 1 int 

Since the Fourier drift scales linearly with acceleration, 
the reduction in sensitivity also scales linearly with ac- 
celeration. For our integration time of 420 s, the critical 
acceleration is a C rit/P ~ 1.7 (where a cr it is in units of m/s 2 
and P is the pulsar period in ms). Our weakest source is 
several times brighter than our detection limit (Figure ^) , 
so a reduction in sensitivity by a factor of several should 
still maintain detectability to pulsations. Thus, a more 
appropriate critical acceleration for the weakest source in 
our list is a cr ;t/P ~ 12 (for P > 1 ms) and a crlt /P ~ 7 
(for P ~ 1 ms). 

Orbits containing millisecond or sub-millisecond pulsars 
that have been spun up from mass transfer from a low- 
mass donor would be expected to be circular. Of the 40 
known pulsars in circular (e < 0.01) binary orbits, the 
largest projected mean acceleration to our line of sight (as- 
suming i = 60°) is Omean/P ~ 3.6 for PSR J1808-3658, an 
X-ray millisecond pulsar with a 2.5 ms period in a 2 hr bi- 
nary orbit around a 0.05 Mq companion (Chakrabarty & 
Morgan 1998). This acceleration is well below the critical 
acceleration a cr it/P even for our weakest source. 

The mean flux density of the sources on our list, how- 
ever, is much higher (Sgio ~ 400 mJy, assuming a = 1.6) 
than our weakest source, which raises the critical accel- 
eration for our typical source. Only very large accelera- 
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tions (a moan /P ^ 85 for P > 1 ms and a moan /P ii 40 for 
P $ 1 ms) would prevent detection of pulsations for our 
typical source. A system such as PSR J1808— 3658 with 
5*610 ~ 400 mJy would still be detectable if it had P ~ 0.3 
ms or if it had an orbital period Pb ~ 15 minutes (but 
not both). Thus it is unlikely that binary motion in a 
sub-millisecond or millisecond pulsar system would be a 
significant source of non-detections for most of our sources. 

We have also estimated the likelihood of serendipitously 
detecting a pulsar not associated with these sources using 
the observed surface density of both normal and millisec- 
ond pulsars in the Galactic plane (Lyne et al. 1998). With 
standard assumptions for a spectral power law index and 
luminosity distribution (Lorimer et al. 1993), we find that 
it is unlikely (< 2% probability) that we would detect any 
pulsars from the chance placement of our 92 beams. 

5. CONCLUSIONS 

No pulsations were detected at 610 MHz from the 92 
polarized, point-like sources that we searched from the 
FIRST and NVSS radio surveys. Sensitivity to sub- 
millisecond pulsations was maintained for DMs less than 
about 500 pc cm~ 3 (without scattering effects), which en- 
compasses the expected DM range for all of these sources. 
We find that several effects which could prevent detection 



(brightness, dispersion smearing, scattering, and beaming) 
are not significant factors here. Scintillation is expected to 
account for only a few of our non-detections and therefore 
cannot be the cause of the majority of our non-detections. 
For a source with a typical flux density in our list, Dopplcr 
motion in a tight binary system would only prevent detec- 
tion if the mean projected line-of-sight acceleration of the 
pulsar were at least an order of magnitude higher than 
those observed in the known population of circular binary 
pulsar systems. We conclude that as a population, these 
sources are unlikely to be pulsars. Given that ~ 10% of 
extragalactic sources in the NVSS survey were identified 
by Han & Tian (1999) as being at least 5% linearly po- 
larized, it is possible that most of our target sources are 
unidentified extragalactic objects. However, the nature of 
these sources is still not certain. 
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A Search for Pulsations in FIRST and NVSS Sources 



Table 1 — Continued 



Source 


FIRST 


a 


8 


Sl400 a 


% polarization b 




Source? 


(J2000) 


(J2000) 


(mjy) 




1150+3020 


y 


11:50:43.890 


+30:20:17.66 


31 


13 


1201+2550 


y 


12:01:25.648 


+25:50:04.55 


20 


24 


1201+3129 


y 


12:01:44.264 


+31:29:03.22 


87 


6 


1220+3111 


y 


12:20:04.656 


+31:11:45.04 


28 


8 


1234+2917 


y 


12:34:54.323 


+29:17:43.93 


434 


9 


1236+3706 


y 


12:36:50.831 


+37:06:02.01 


62 


9 


1242+2721 


y 


12:42:19.687 


+27:21:57.32 


69 


10 


1251+3643 


y 


12:51:24.132 


+36:43:57.27 


29 


8 


1334+3434 


y 


13:34:26.833 


+34:34:25.11 


48 


7 


1343+2903 


y 


13:43:24.006 


+29:03:57.55 


22 


17 


1414+4022 


y 


14:14:40.585 


+40:22:25.75 


43 


6 


1426+4035 


y 


14:26:58.101 


+40:35:38.36 


28 


9 


1434+3805 


y 


14:34:46.988 


+38:05:14.87 


149 


9 


1458+3720 


y 


14:58:44.704 


+37:20:22.15 


210 


6 


1508+2818 


y 


15:08:08.312 


+28:18:13.51 


76 
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1547+3954 


y 


15:47:40.147 


+39:54:38.48 


128 


15 


1606+2709 


y 


16:06:16.249 


+27:09:28.69 


29 


7 


1609+2628 


y 


16:09:50.978 


+26:28:38.72 


17 


15 


1618+2931 


y 


16:18:27.685 


+29:31:17.99 


29 


8 


1635+3751 


y 


16:35:53.071 


+37:51:54.59 


47 


6 


2321-1758 


n 


23:21:02.411 


-17:58:22.09 


17 
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a NVSS catalog total intensity 1400 MHz peak flux density. 

b NVSS catalog percent linear polarization from 1400 MHz peak flux densities. 



Crawford ct al. 
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Table 2 
Observed Test Pulsars 



Name 


P 


DM 


5l400 a 


Seoo b 


S/N 




(ms) 


(pc cm -3 ) 


(mjy) 


(mjy) 




PSR B0329+54 


714.52 


26.8 


203 


785 


725 


PSR B1937+21 


1.56 


71.0 


16 


100 


46 


PSR J2145-0750 


16.05 


9.0 


10 


30 


36 



a Catalog 1400 MHz flux density (Taylor, Manchester, & Lyne 1993). 
b Catalog 600 MHz flux density (Taylor, Manchester, & Lyne 1993). 
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A Search for Pulsations in FIRST and NVSS Sources 
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Fig. 1. — Pulsar sensitivity curves for DMs of 0, 100, and 500 pc cm~ 3 , assuming a 5% intrinsic pulsed duty cycle. The dashed horizontal 
line at 58 mjy is the flux density of our weakest source at 600 MHz, assuming a typical pulsar spectral index of a = 1.6. All of our sources 
are expected to have DM < 100 pc cm~ 3 (indicated by the dashed-dotted line). 



